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Open dielectric resonators of different shapes are widely used for the manufacture of microlasers. 
A precise determination of their resonance frequencies and widths is crucial for their design. Most 
microlasers have a flat cylindrical geometry, and a two-dimensional approximation, the so-called 
method of the effective index of refraction, is commonly employed for numerical calculations. Our 
aim has been an experimental test of the precision and applicability of a model based on this ap- 
proximation. We performed very thorough and accurate measurements of the resonance frequencies 
and widths of two passive circular dielectric microwave resonators and found significant deviations 
from the model predictions. From this we conclude that the model generally fails in the quantitative 
description of three-dimensional dielectric resonators. 

PACS numbers: 42.55.Sa, 05.45.Mt, 42.60.Da 



I. INTRODUCTION 

Open dielectric resonators are used in a large vari- 
ety of applications, ranging from radio frequency and 
millimeter- wave applications [1, 2j to microlasers ^[4l|5]. 
Therefore, accurate model predictions for their spectra 
and field distributions are of great interest. Especially 
microlasers have received much attention lately in opti- 
cal telecommunication, as sensors or as billiard models 
[2j |5l [H [7] . These devices typically consist of a flat cylin- 
drical dielectric resonator with a cross section of arbitrary 
shape which contains the active medium. The resonators 
are usually made of semiconductor [H |8l |9] or organic ma- 
terials [10l[TTl[T2 and are sandwiched between two media 
of lower index of refraction like air or a substrate. The 
exact shape of the resonator is important in view of the 
applications, because it determines the emission proper- 
ties like the directionality of radiation and the quality 
factors of the resonances [131 El IE] • 
In general, even simple geometries like a flat dielectric 
disk with a height much smaller than the planar ex- 
tension cannot be solved analytically. Since the nu- 
merical solution of the three-dimensional (3D) vectorial 
Maxwell equations describing such resonators is compli- 
cated and computationally demanding, suitable approx- 
imations are favorable. One widely used approximation 
is the reduction of the full 3D Maxwell equations to a 
two-dimensional (2D), scalar Helmholtz equation by in- 
troducing a so-called effective index of refraction neff (see 
e.g. [HI [171 [18^ 19 ). This 2D-approximation seems natu- 
ral due to the fiat shape with large extension in the plane 
of the microlasers, but to our knowledge its validity and 
precision has never been rigorously tested. The aim of 
the work presented here is thus the comparison of the ex- 
perimentally measured resonance frequencies and widths 
with those calculated using this 2D-approximation, called 



the neff-model in the following. 

Although motivated by microlasers working in the in- 
frared to the optical spectrum, our experimental setup 
consists of a fiat cylindrical microwave resonator made of 
Tefion (from the company Griinberg Kunststoffe GmbH) . 
Microwave experiments have distinct experimental ad- 
vantages, and therefore are commonly used for the inves- 
tigation of 2D quantum billiards [20l [21] and also 2D di- 
electric resonators [22 . The results from the microwave 
resonators can, however, be directly applied to micro- 
cavities by scaling. The plan of the present article is the 
following. In section [n| the concept of an effective index 



of refraction is introduced and section III explains how 
dielectric resonators are modeled with it. The experimen- 
tal setup is detailed in section [TVl and the experimental 
data are compared to the model calculations for two disks 
of different thickness in sectio ns [V| and |VI[ Finally, the 
results are discussed in section IVIII 



II. EFFECTIVE INDEX OF REFRACTION 

The basic idea for the 2D-approximation is to treat the 
bulk of the resonator as a dielectric slab waveguide [23] . 
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FIG. 1: Geometry and notations for the infinite dielectric 
slab waveguide. The dielectric slab with index of refraction 
n and thickness h is extended indefinitely in the x-y-plane 
and surrounded by media with index of refraction ni and n2, 
respectively. 
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FIG. 2: Effective index of refraction rieff with respect to kb 
for n = 1.434. The sohd hnes are the TE-modes, the dashed 
hnes the TM-modes of various z-excitations. The dotted hnes 
are the indices of refraction of Teflon and air, respectively. 



In this section we will consider an infinite dielectric slab 
waveguide, that is an infinite dielectric plate of thickness 
b with index of refraction n surrounded by media with 
lower indices of refraction ni^2 (see Fig. [l]), and treat the 
cylindrical sidewalls of the resonator in the next section. 
The surrounding media are assumed to be air (ni,2 = 
1) in the following. The waveguide is described by the 
vectorial Helmholtz equation 



E 
B 



= 



(1) 



where k = co/c is the vacuum- wavenumber, co the angu- 
lar frequency and c the speed of light. The Helmholtz 
equation ([T]) can be simplified by separation of the vari- 
able parallel to the cylinder axis, 2:, from the x and y 
variables. Different modes in the slab waveguide can be 
distinguished by their polarization and their excitation 
perpendicular to the plane of the slab. Due to the slab 
geometry all field modes are either related to (TM- 
polarization) or to Bz (TE-polarization) [24 . The gen- 
eral ansatz for E^ respectively is 
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and the constants are determined from the boundary 
conditions. The wave function ^ satisfies the scalar, two- 
dimensional Helmholtz equation 



(A + 7^)^ = 0, 



(3) 



with 7 being the horizontal component of the wavevector. 
The vacuum- wave numbers k and its vertical components 
kz inside the dielectric medium and Qz outside are related 
by the dispersion relation 



k' 



(4) 



It should be noted that qz is real for all solutions which 
correspond to modes confined inside the dielectric slab 
by total internal reflection with evanescent fields outside. 



The continuity condition for in?{r)Ez and at the in- 



terfaces at z = ±6/2 for TM-modes and Bz and 
TE-modes yields the condition [16] 
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Expressing the horizontal component 7 of the wavevector 
in terms of the effective index of refraction neff, defined 



as 



neE = l/k, (6) 
in Eq. Q and Eq. ([5| leads to the dispersion relation 



kb/2 : 



arctan v 
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for TM^ 
for TEq 



(7) 



for neff with = 0, 1, 2, . . . denoting the order of exci- 
tation in z-direction. Since in the following only TMq 
and TEo modes are considered, we skip the index. Obvi- 
ously, neff depends only on the index of refraction n and 
on kb oc 6/A, that means on the ratio of the slab's thick- 
ness b to the wavelength A. A plot of neff with respect to 
kb is shown in Fig. |2] The effective index of refraction is 
always 1 < neff < n. Modes with z-excitation ( emerge 
at the cutoff kb ■ 
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III. 



APPLICATION TO DIELECTRIC 
RESONATORS 



The solution of neff is now inserted into Eq. ([3|, lead- 
ing to the two-dimensional scalar Helmholtz equation 
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k^)^ = 



(8) 



with ^ corresponding to Ez {Bz) for TM- (TE-) modes. 
This equation correctly describes the propagation of elec- 
tromagnetic waves inside an infinite dielectric slab. The 
next step is to incorporate the cylindrical sidewalls of 
the dielectric resonator in a plane perpendicular to the 
cylinder axis. In the 2D approximation this is achieved 
by considering the vertical walls as a part of an infinite 
dielectric cylinder with n = neff and imposing the corre- 
sponding boundary conditions [24], i.e. that fields inside 
and outside the resonator obey the Helmholtz equation 



in, out 



"eff 
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Here S is the domain of the resonator in a plane perpen- 
dicular to the cylinder axis, and the conditions along the 
boundary dS of S 
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are imposed, with /i equal to 1 (l/nl^) for TM- (TE-) 
modes and n being the unit vector normal to the sur- 
face. The set of equations (9| and (10) constitutes the 
quintessence of the neff-model. It correctly describes a 
two-dimensional dielectric resonator, i.e. the electromag- 
netic field is homogeneous in the z-direction, with in- 
dex of refraction equal to neff, however, not a fiat but 
three-dimensional resonator. Indeed the full 3D Maxwell 
equations lead to additional boundary conditions if the 
electromagnetic fields depend on z [T9^, ^ , which couple 
the TM- and TE-polarizations. Analytical calculations 
incorporating these have to our knowledge not yet been 
performed. As has already been stated above, the aim 
of the present work is to test experimentally the appli- 
cability of the neff-model and to understand the order of 
magnitude of deviations from the experiment. 
The solutions of Eqs. ([9| and ([lO| for a circular resonator 
with radius R are given in cylindrical coordinates r and 
cp as 



suspensions 



^in(r,(/:)) 



,(neff/cr)e 



(11) 



inside the disk and as 

*out(r,^) = *i°lHW(fcr)e±*'"^ (12) 

outside. Here, Jm(^) is a Bessel- function of the first kind, 
H^^(x) a Hankel-function of the first kind and m the 
azimuthal quantum number. All modes with m > are 
doubly degenerate. Applying the boundary conditions 
(Eq. (10)) leads to the quantization condition [26 



JmiriefikR) H^^(A;i?) 
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where J'm(^) and H'^^(x) are the derivatives with re- 
spect to X. For each azimuthal quantum number m, 
there is an infinite series of complex solutions km^ur'^ 
nr = 1, 2, . . . is the radial quantum number. Since a di- 
electric resonator described by equations (l9|) and ( 10 ) is 



an open system, its modes have losses due to radiation. 
These are determined by the imaginary part of km,nr^ 
and the quality factor Q of a mode is defined as 



2lm{km,nr.)' 



(14) 



Because of these losses, the solutions of Eq. ([9| are called 
quasi-bound modes. They are identified as the poles of 
the scattering matrix S describing the measurement pro- 
cess [3 — microwave power is coupled into the resonator 
via an antenna, thereby exciting modes, and coupled out 
via the same or another antenna — in the form of reso- 
nances in the frequency spectrum. 

The quantization condition Eq. ([l3| is solved numeri- 
cally, taking into account the dependence of neff on (the 
real part of) k. 
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FIG. 3: Schematic side view of the experimental setup. The 
Teflon disk is hanging on three metal suspensions. Two dipole 
antennas protruding from an RF cable are placed close to the 
rim of the disk on opposite sides. 





FIG. 4: Sketches of the two antenna types. A dipole an- 
tenna is shown on the left, a curved antenna on the right. 
The cylindrical sidewall of the Teflon disk is shown in gray in 
the background. Both types of antennas are placed directly 
alongside the sidewall of the disk to obtain good coupling to 
the resonator. 



IV. EXPERIMENTAL TECHNIQUE 

The experimental setup is sketched in Fig. [3] The 
Tefion disk is hanging on three metal suspensions at- 
tached to the disk at the corners of an equilateral trian- 
gle, but otherwise surrounded only by air. Two antennas 
are put on opposite sides of the disk. The whole setup 
is placed in a thermostat to keep it at a fixed temper- 
ature. It should be noted that the perturbation of the 
resonator due to the attached suspensions leads to neg- 
ligible changes of the resonance frequencies by less than 
2 % of the mean resonance spacing. The perturbations 
due to the bending of the resonator under its own weight 
and the antennas are of the same order of magnitude. 
Two different types of antennas (see Fig. [4]) were used: 
vertical dipole antennas, which excite mainly TM- modes, 
and so-called curved antennas, which excite mainly TE- 
modes. The antennas may slightly lift the degeneracy of 
the modes. A vectorial network analyzer (PNA N5230A 
by Agilent Technologies) is used to measure the scatter- 
ing matrix element 6'2i(/), where the modulus squared 
of S2i{f) is the ratio 



\S2l{f)\' 



(15) 



between the power Pin coupled in by antenna 1 and the 
power Pout coupled out by antenna 2 for a given fre- 
quency / = a;/27r. 
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FIG. 5: Frequency spectrum of disk A measured with dipole antennas. The transmission amphtude, i.e. the modulus |5'i2|, is 
shown with respect to the frequency. In the magnified part of the spectrum, the resonances are labeled with TM/TE (m,nr) 
to indicate their polarization as well as their azimuthal and radial quantum numbers m and rir, respectively. Two series of 
resonances can be seen here: The broader and larger resonances correspond to modes with TM-polarization and radial quantum 
number rir = 1, the sharper and smaller resonances to modes with TE-polarization and Ur — 1. Resonances with rir > 1 can 
also be seen at higher frequencies. 




7.94 GHz: TM (50, 1) 



13.80 GHz: TM (92, 2) 
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FIG. 6: (Color online) Measured intensity distributions of three TM-modes with quantum numbers (m,nr). A mode with 
azimuthal quantum number m and radial quantum number rir has 2m maxima in azimuthal direction and rir rings. Shown are 
modes of whispering gallery type, as are all other identifiable modes. Therefore, intensities were measured only in the outer 
ring of the resonator. 



Two circular disks of different thickness h made of Teflon 
were used in the experiments to investigate the depen- 
dence of the resonance frequencies on the thickness. 
Disk A has a radius of = 274.8 mm and a thickness 
of 6 = 16.7 mm, and disk B of = 274.9 mm and of 
h = 5.0 mm. A frequency of 10 GHz corresponds to 
kR ^ 57.6, and kb ^ 3.5 (disk A) and kb ^ 1.0 (disk B), 
respectively. The index of refraction was measured using 
a split-cylinder resonator [27, 28 and is n = 1.434 ± 0.01 
for disk A and n = 1.439 ± 0.01 for disk B. An example 



of a measured frequency spectrum of disk A (with dipole 
antennas) is shown in Fig. [5] It displays a superposition 
of several series of almost equidistant resonances, each 
corresponding to modes with a fixed polarization and ra- 
dial quantum number rir and ascending azimuthal quan- 
tum number m. The resonance spacing for each of these 
subspectra is typically 120-130 MHz. This is illustrated 
in the lower part of Fig. [5j Within each subspectrum, 
the width of the resonances decreases with increasing az- 
imuthal quantum number. This can be explained within 
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the ray picture [26]: A higher azimuthal quantum num- 
ber corresponds to a higher angular momentum, and this 
to rays with a larger angle of incidence at the boundary, 
implying lower radiation losses. As a consequence, differ- 
ent subspectra are distinguishable only above a certain 
frequency and contribute below only to the background. 
The dependence of the amplitude of a resonance on the 
type of antenna used gives a hint at the polarization of 
the mode. It was determined in addition with a pertur- 
bation technique: A metal plate was introduced parallel 
to the Teflon disk (with varying distance) leading to a 
shift of the resonance frequencies. Due to the different 
boundary conditions for TE- (Dirichlet for B^) and TM- 
(Neumann for E^) modes at the metal plate the former 
are shifted to lower, the latter to higher frequencies with 
decreasing distance between the metal plate and the disk. 
The quantum numbers were determined from the inten- 
sity distributions, which were measured with the pertur- 
bation body method [29]. A cylinder made of magnetic 
rubber [30] was used as perturbation body and moved 
along the surface of the disk, its height of 8 mm and 
diameter of 4 mm being small compared to a vacuum- 
wavelength of 30 mm at 10 GHz. Then the positioning 
of the perturbation body on the disk leads to a shift of 
the resonance frequency which is proportional to the elec- 
tric field intensity at its position. Three examples are 
shown in Fig. [6| All the measured modes are of whis- 
pering gallery mode type. As a result the resonances are 
only slightly perturbed by the suspensions, since they are 
located well inside the caustic of the whispering gallery 
modes. Accordingly, the intensities were measured only 
in the outer part of the disk. With the knowledge of the 
polarization and the quantum numbers, the measured 
resonance frequencies can be compared with those com- 
puted based on the neff-model (solutions of Eq. (13)). 



V. COMPARISON OF MODEL AND 
EXPERIMENT FOR DISK A 

Figure [T] shows the difference between the measured re- 
sonance frequencies /expt of disk A {b = 16.7 mm thick) 
and those calculated with Eq. ( [l3| ), /caic- In Fig- |7(a)| 
(TE- modes), there are two series of data points for each 
radial quantum number because the degeneracy of 
the modes is lifted by the curved antennas used in the 
measurement. The scattering of the data points of about 
±5 MHz in both graphs stems from problems with the de- 
termination of the resonance frequencies, either because 
the resonances are badly shaped (at lower frequencies) 
or because of overlapping resonances (at higher frequen- 
cies). Some resonances are not detectable due to the 
overlap with others. The deviation between the measured 
and computed resonance frequencies is generally less than 
1 %. Still, these deviations of the model calculations from 
the experimental data must be considered significant: A 
deviation of 60 MHz is about half the resonance spac- 
ing between resonances with the same n^, and thus it is 
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FIG. 7: Difference between measured (/expt) and calculated 
(/caic) resonance frequencies with respect to /expt for disk 
A. The different symbols correspond to the different radial 
quantum numbers (x: rzr = 1, o: rir = 2, +: rir — 3). 
(a) TE-modes: The range of azimuthal quantum numbers for 
rir 1 is m = 37-148. For each rir, there are two series of 
data points due to the break-up of the degenerate modes by 
the curved antenna used in the measurement. The frequencies 



of the unperturbed system are approximately in between, (b) 
TM-modes: The range of azimuthal quantum numbers for 
rir = 1 is m = 39-150. The measurement was done with 
dipole antennas, and no break-up of degenerate modes was 
observed. 



impossible to correctly identify the different resonances 
just by a comparison with the model calculations. The 
difference between the calculated and the measured re- 
sonance frequencies decreases with increasing frequency 
for both polarizations and seems to reach a finite value 
in both cases. Interestingly, this value is different for 
TE- and TM-modes. Furthermore, the magnitude of the 
deviations slightly depends on the radial quantum num- 
ber, especially for the TE-modes. All this hints at a 
systematic failure of the model, even in the semiclassical, 
i.e. high frequency limit, although the decrease of the 
deviations with increasing frequency indicates that the 
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FIG. 8: Difference between measured (/expt) and calculated 
(/caic) resonance frequencies with respect to /expt for the TE- 
modes of disk A. Results from three different calculations with 
three different indices of refraction n — 1.432 (o), 1.434 (x) 
and 1.436 (+) are shown only for resonances with rir — 1 
for the sake of clarity. The difference /expt — /caic depends 
strongly on n in the semiclassic limit, i.e. for high frequencies, 
but only very slightly for low frequencies. 



neff-model is more accurate in the semiclassical regime. 
It should be noted that the magnitude of deviations be- 
tween model and experiment is extremely sensitive to the 
exact value of the index of refraction, which has been de- 
termined with an uncertainty An = ±0.01. In order to 
illustrate that this accuracy is insufficient for an exact 
determination of the deviations, /expt — /caic is plotted in 
Fig. |8]for calculations with three values of n, which dif- 
fer by even less than An. Only data points with n^ = 1 
are shown, the data points for resonances with n^ > 1 
go along with them. The three calculations have a very 
different behavior in the semiclassical limit, but roughly 
agree for low frequencies. 

In a next step we considered the index of refraction n as a 
fit parameter — let's call it n — in Eq. (13), where n en- 



ters implicitly via neff, and varied it such that Eq. (13) 
yields the measured resonance frequencies. If the neff- 
model provides a good description, the resulting values 
of n should scatter around the actual value of n. The 
results are shown in Fig. [9) Almost all data points (ex- 
cept the ones at low frequencies) lie inside the error band 
n — An, but they form three distinct curves for the three 
values of the radial quantum number n^. This provides 
further evidence for the failure of the neff-model: If the 
model were correct, the values of n would form a single 
line along the actual index of refraction n of the disk for 
both polarizations and all three values of n^. It is known 
from literature (and was confirmed experimentally) that 
Teflon has negligible dispersion in this frequency range. 
In conclusion the observed deviations between the model 
and the experiment cannot be attributed to a badly de- 
termined index of refraction, as it is impossible to achieve 
agreement between /expt and /caic in the whole frequency 
range by choosing a fixed value of n. This is also true 
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FIG. 9: Values of the index of refraction ft for which Eq. ( 13 ) 



yields the measured resonance frequencies /expt with respect 
to the measured resonance frequency for disk A. The upper 
graph (a) shows the result for the TE-modes, the lower one 



(b) that for the TM-modes. The different symbols denote the 
different radial quantum numbers (x:nr = l, o:nr = 2, +: 
Ur — 3). The solid line is the real index of refraction n of 
the disk, the dashed line n — An. The systematic deviation 
of the data points from the measured n shows the failure of 
the rieff-model. 



for the other two parameters, the radius R and the thick- 
ness 6, or combinations thereof. Thus we can exclude the 
possibility of badly determined parameters and inaccura- 
cies in the measurement of the resonance frequencies and 
state our main result: the neff-model does not correctly 
describe the measured resonance frequencies. 
Finally, the measured and calculated resonance widths 
(FWHM) are compared. The experimental resonance 
widths Texpt are obtained by fitting Lorentzians to the 
measured spectrum. They consist of three terms. 



Texpt — r^acj -|- rai3s -|- Tg 



(16) 



where Frad describes the losses due to radiation, Fabs 
the losses due to absorption in the Teflon and Fant the 
loss of power due to the coupling to the antennas. The 



7 




7 8 9 10 11 12 13 

/expt (GHz) 

FIG. 10: Measured (Fexpt) and calculated (Fcaic) resonance 
widths for TM-modes of disk A. Both Fexpt and Fcaic are plot- 
ted as function of the measured resonance frequency /expt- 
The different symbols denote the measured widths for differ- 
ent radial quantum numbers (+: rir = 1, x: rir = 2). The 
calculated widths are plotted as curves (solid line for rir = 1, 
dashed for rir = 2) instead of data points to guide the eye. 



calculated resonance widths 



Fcaic = -2cIm(^)/(2^) 



(17) 



include only the losses due to radiation (Frad)- The mea- 
sured and calculated widths of the TM-modes for disk A 
For low frequencies (up to 8 GHz 
2), the calculated 
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are shown in Fig. 
for rir = 1 and up to 10.5 GHz for 
widths are up to twice as large as the measured widths. 
Since Fcaic does not account for absorption and antenna 
losses, the actual difference is even larger. For higher fre- 
quencies, the measured widths are larger than the calcu- 
lated ones and saturate at a value of about 4 MHz. This 
saturation is due to absorption in the Teflon material and 
coupling to the antennas, and is approximately indepen- 
dent of the frequency. A precise comparison with the 
calculated widths is generally not possible because the 
radiation losses F^ad cannot be extracted from the mea- 
sured widths Fexpt- Nonetheless it can be stated that the 
widths Fcaic predicted by the neff-model are too large at 
least in some frequency ranges within the range of accu- 
racy of the index of refraction n. It should also be noted 
that Fcaic does not depend as sensitively on n as the 
resonance frequencies. For TE-modes, the same general 
trend for Fgxpt and Fcaic is found, although the difference 
between Fexpt and Fcaic is not as pronounced as for the 
TM-modes. 



VI. COMPARISON OF MODEL AND 
EXPERIMENT FOR DISK B 

Next, we will compare the data of the second, thin- 
ner disk B (6 = 5.0 mm) with the calculations based on 
the neff-model. The discrepancy between experiment and 
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FIG. 11: Difference between measured (/expt) and calculated 
(/caic) resonance frequencies with respect to /expt for disk 
B. The different symbols correspond to the different radial 
quantum numbers (x: rir = 1, o: rir = 2, +: rir = 3). 
The T E-m odes measured with curved antennas are plotted in 
graph |(a)| and have a range of azimuthal quantum numbers 
m — 64-188 for resonances with rir — 1. The TM-modes 



shown in graph (b) were measured with dipole antennas and 
have azimuthal quantum numbers m = 97-204 for rir = 1. 
The frequency range of identifiable TE- and TM-modes differs 
due to the different quality factors and types of antennas used. 



model becomes apparent when comparing measured and 
calculated resonance frequencies in Fig. [TT] (like in Fig. [7| 
for disk A). The differences /expt — /caic seem to be larger 
for disk B, but are of the same order of magnitude. How- 
ever, in contrast to those for disk A, they increase with 
increasing excitation frequency. In fact, this behavior de- 
pends sensitively on n; changing the value of n within the 
range of accuracy leads to differences increasing, decreas- 
ing or reaching a finite value with increasing /expt (see 
also Fig. [8]). In any case the deviations are larger than 
for disk A at least for frequencies up to 20 GHz. 
Like in Fig. |9j the index of refraction h needed to repro- 
duce the measured resonance frequencies with Eq. (13) 
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FIG. 12: Index of refraction n required to reproduce the mea- 
sured resonance frequencies with Eq. (13) as a function of 
the resonance frequency /expt for disk B. Each symbol corre- 
sponds to a radial quantum number (x : rir = 1, o: nr = 2, 
+: rir = 3), the solid line denotes the index of refraction 
n = 1.439 of the disk and the dashe d lin e n — An the range of 
accuracy of its determinati on. P art | (a) [shows the data points 
for the TE-modes and part (b) those for the TM-modes. The 



trend of the curves is comparable to Fig. [9] but the deviations 
from n are even larger. 



is shown for = 1,2,3 in Fig. [12] with respect to the 
resonance frequency /expt- Qualitatively, the curves are 
similar to those in Fig. [9j n grows with increasing fre- 
quency and seems to converge to some value depending 
on the polarization, and h has a weak dependence on the 
radial quantum number of the modes. Compared to 
the results for disk A shown in Fig. [9] the deviation of n 
from n seems to be even larger for the thinner disk B, and 
those between the curves for different are more pro- 
nounced. Again, the systematic failure of the neff-model 
is clearly visible, and the discrepancies seem to be even 
larger. 

The comparison of calculated and measured resonance 
widths for disk B is plotted in Fig. [13] Like in Fig. [Toj 
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FIG. 13: Measured (Fexpt) and calculated (Fc 
widths for the TE-modes of disk B. Both F, 



expt 



20 



resonance 
and Fcaic 



are plotted with respect to the measured resonance frequency 
/expt- The different symbols denote the measured widths for 
different radial quantum numbers (x: rir = 1, o: rir = 2, +: 
rir = 3). The calculated widths are plotted as curves (solid 
line for rir = 1, dashed for rir = 2 and dot-dashed for rir = 3) 
instead of data points to guide the eye. 



the measured widths Fexpt are significantly smaller than 
the calculated widths Fcaic for the resonances with lower 
azimuthal quantum numbers, although the difference is 
not quite as large as for the case of disk A. With increas- 
ing frequency, the measured widths saturate at a value 
of about 7 MHz, and a comparison with the calculated 
widths is not possible because of the additional losses due 
to absorption and the antennas. The difference between 
Fexpt and Fcaic IS Smaller in the case of TM-modes, and 
gets smaller for both cases if a larger index of refraction 
n is assumed in the calculations. Still, it is definite that 
the calculated widths are too large at least for some re- 
sonances. 



VII. CONCLUSIONS 

We have measured the resonance frequencies and 
widths of two different flat cylindrical dielectric mi- 
crowave resonators. The quantum numbers and the po- 
larization of the corresponding modes were identified. 
These data were used to test the neff-model, a two- 
dimensional approximation for flat dielectric resonators 
with large planar extension normally used for the model- 
ing of e.g. microlasers. A microwave resonator was chosen 
as testbed for several reasons: It is easy to handle due to 
its macroscopic dimensions which are known with high 
precision, the measurement of intensity distributions en- 
ables a detailed understanding of the spectrum and the 
identification of individual resonances, and it is a passive 
system which means that additional shifts of resonance 
frequencies and widths due to the lasing process in an 
active medium [31, ^32^ can be excluded. 
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It was shown that the rieff-model fails to correctly pre- 
dict the resonance frequencies in a systematic way, and 
it has been checked carefully that the deviations between 
experiment and calculations are not due experimental in- 
accuracies caused by e.g. the antennas or the suspensions. 
Also, the resonance widths are clearly overestimated at 
least in the lower frequency range. This is important be- 
cause the widths (of the passive cavity) determine the las- 
ing threshold in microlasers and play a crucial role with 
regard to mode-competition [STJ [32] . Although the calcu- 
lations based on the neff-model yield at least the correct 
order of magnitude for both resonance frequencies and 
widths, a detailed understanding of the spectrum with 
help of these calculations is impossible. Furthermore, the 
accuracy of the model is not under control: The magni- 
tude of the deviations between model and experiment 
cannot be precisely determined due to the uncertainty in 
the index of refraction, and depends on the dimensions 
of the cavity and the frequency in a non-trivial way. It 
seems that deviations are larger for smaller values of kb^ 
i.e. thinner resonators. We believe that the main reason 
for the failure of the neff-model is the incorrect formu- 
lation of the boundary conditions (as in Eq. (10)). In 



the rigorous treatment of the 3D Maxwell equations the 
z-components of the electric and magnetic fields are con- 
nected by the boundary conditions and cannot be con- 
sidered separately as done in the neff-model [19l [25] . 
The large deviations between model and experiment ob- 
served for the circular disk are not due to the localization 
of the (whispering gallery) modes close to the boundary, 
and not expected to be smaller if modes supported by 
the whole area of the disk were considered. In fact the 
modes with lower n^ are localized closer to the bound- 
ary (see Fig. |6]), but the deviations are smaller than for 
modes with higher n^. Even though the circular res- 
onator is a special case and the test of the neff-model with 



a non-circular geometry remains an open problem, it is 
still worthwhile to consider it: the quasi-bound states of 
many other non-circular resonators are also of whispering 
gallery mode type. Moreover, the circular resonator has 
many important applications by itself 2^ • One of these 
is the precise measurement of the index of refraction over 
a wide range of frequencies which, however, relies on a 
rigorous model for the resonance frequencies. 
In summary, the comparison of the measured and com- 
puted resonance frequencies and widths clearly attests 
the need of an improved neff-model, which takes into ac- 
count the boundary conditions at the cylindrical side- 
walls. We hope that this work will stimulate further re- 
search in this direction to obtain a reliable 2D model for 
the computation of the resonance frequencies and widths 
of fiat dielectric resonators. Although the numerical so- 
lution of the full three-dimensional Maxwell's equations 
is feasible for the circular resonator [33*, '3T, "35 , it is com- 
putationally demanding, and even more so for other (viz. 
chaotic) geometries. Still, 3D numerical calculations for 
the circular resonator would be helpful to validate our 
data and investigate the accuracy of the neff-model in a 
broader range of aspect ratios b/R and indices of refrac- 
tion n. 
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